P. Filip 



Elliptic Flow in Collisions of Deformed Nuclei 



Nuclear geometry and elliptic flow in 
collisions of deformed nuclei 



P. Filip 

Institute of Physics, Slovak Academy of Sciences, Bratislava 

... I Introduction 

o\ ■ 

Elliptic flow measurement jU [2] has become a very important method for determining prop- 
erties of dense strongly interacting matter created in ultra-relativistic heavy ion collisions exper- 
iments. Being theoretically predicted in 1992 [3] the elliptic flow has been observed in heavy ion 
h- » i collisions at AGS, SPS and RHIC % 

The elliptic flow is generated during the expansion of compressed QCD matter into vacuum 
from the initial asymmetrical volume created in non-central collisions of nuclei. During the 
process of expansion (which can be described by hydro dynamical [3] or rescattering [3] models) 
an asymmetry in the azimuthal distribution of transverse momenta of particles is generated. 

Experimentally measured azimuthal distribution of particles in non-central heavy-ion collisions 
exhibits second-order oscillation of the form = a[l + 2t>2 • cos(20)], where parameter v<i 
i characterizes the strength of the elliptic flow effect [lj. Quantitative comparisons of the initial 
spatial eccentricity e and final elliptic flow strength w 2 allow one to make implications on the 
equation of state of strongly-interacting QCD matter [5]. 

In this contribution we investigate the elliptic flow behaviour in relativistic collisions of de- 
^ | formed nuclei. Using Optical Glauber Model simulation |7j we predict large fluctuations of the 
; elliptic flow (at given collision centrality) in collisions of deformed nuclei. We also suggest that 
behaviour of the elliptic flow in the most central collisions of deformed nuclei is different for nuclei 
^ . with oblate and prolate deformations. 

O 



II Elliptic flow in collision of spherical nuclei 



c5 ' Initial eccentricity e of the compressed QCD matter in collisions of spherical nuclei (e.g. Pb 207 ) 
is simply related to the size of impact parameter b (a distance between centers of colliding nuclei 
in transversal plane). Collisions of spherical Pb nuclei with fixed impact parameter (e.g. 6 = 3 fm) 
will exhibit elliptic flow values fluctuating around some average value (i^). Fluctuations of V2 
values in such collisions with fixed impact parameter originate e.g. from the fluctuations of initial 
participant eccentricity which occur due to varying positions of individual nucleons [8J inside the 
interaction volume. Thus even in the most central collisions of spherical nuclei (b = Ofm) the 
eccentricity e fluctuates and on average (e) > 0. 

These effects can be studied in Monte Carlo Glauber Model simulations [9] where individual 
positions of nucleons are determined for each collision event. In the optical Glauber model such 
effects are neglected. 
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In the next sections we predict event -by-event fluctuations of the elliptic flow v 2 due to initial 
eccentricity e fluctuations originating from the deformed shape of colliding nuclei. Such v 2 fluc- 
tuations (present only in collisions of deformed nuclei) are generated additionally, on top of v 2 
fluctuations observed in collisions of spherical nuclei. 



Ill Elliptic flow in collisions of deformed nuclei 

In collisions of deformed (oblate or prolate) nuclei the initial eccentricity e of the interaction 
zone strongly depends on the orientation of nuclear spin (axis of the ellipsoid). For example, in 
central (b = Ofm) collision of two prolate nuclei (e.g. Ho 165 ) the eccentricity can be large if spins 
of Ho 165 nuclei are parallel to each other and orthogonal to the beam axis, or negligibly small 
for spins of both Ho 165 nuclei being parallel to the beam axis. Eccentricity e thus depends on 
azimuthal and polar angles <ft, 9 of nuclei colliding: 

e[b] = f(9 1 ,<f> l ,9 2 ,<f> 2 ) (1) 

Using the optical Glauber simulation [7] one can calculate number of participating nucleons 
iVpart, number of nucleon-nucleon collisions N co \\ and eccentricity e for any given value of impact 
parameter b and angles 0i, 9\, 4> 2) 9 2 (see Fig.l). 




Figure 1: Contour-plot of [ e ; -/V par t] values (left) obtained with optical 
Glauber simulation of Ho 165 +Ho 165 collisions using fixed impact parameters 
b = 0, 1, 2, . . . 12 fm, and the shape of Hoi/2_ nucleus with f3 2 = 0.3 (right). 



III.l Prolate nuclei 

Let us consider rare-earth nucleus Ho 165 . This element has a single stable isotope with prolate 
ground-state deformation j3 2 = 0.3 [TU]. Using deformed Woods-Saxon density [UJ the eccentricity 
of the interaction zone 

(2) 



<7 2 4- CT 2 

"j 1 x 
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(here a xy = (x ■ y) — (x)(y)) can be evaluated (e.g. using N co u[x,y] density in transversal plane) 
for all combinations of independent angles 61,62 E (0,90) and (f>i,(f>2 E (0,180). This has been 
done in 15 degree steps for impact parameters b = 0,1, 2,... 14 fm (see Fig.l). Every given 
collision event with angles 61,62 has been weighted by spherical angle factor sin^) sin^) and 
impact parameter factoi0 2n ■ b to account for the probability of the collisions as they occur in 
unpolarized experiments. In the resulting histogram of (s; N part ) values (as shown in Fig.l) one 
observes large fluctuations of eccentricity e for collisions with fixed impact parameter b. Also 
number of participants varies at given impact parameter value due to variation of angles 61, 6 2 
and 0i, 2 . 



III.2 Oblate nuclei 



-Au collisions assume the gold nucleus AuJ 97 , to be 



-21 



3/2H 

2 [12J is predicted 



Almost all simulations of relativistic Au 
spherical. However, this nucleus with quadrupole moment Q ~ 0.55 • 10~" cm 
to have slight oblate deformation in the ground state with 02 ~ —0.13 pU]. In order to investigate 
possible consequences of the predicted oblate deformation of Au 197 nucleus in heavy ion collisions 
we assume it to be deformed with /3 2 = —0.131 [TU]. The result of our optical Glauber model 
simulation is shown in Fig.2. One observes fluctuations of eccentricity and number of participants 
at given fixed impact parameters b = 0, 1,2, ...fm. For spherical nuclei the optical Glauber 
simulation would give a single value of the eccentricity and the number of participants for any 
given impact parameter b. 



"^£eiKrlciTy~pNtJsTfT 



|EccerHiicity [ N_part ] [ 




IN.part] 



Figure 2: Contour-plot (left) of [ e ; -/V part ] values obtained for Au 197 +Au 197 collisions at 
b = 0, 1, 2, ... 12 fm assuming @2 — —0.13. Error bars in the Ecc[iVp art ] plot (right) show 
the width of eccentricity fluctuations at given fixed impact parameter values. 



One can conclude that ground-state deformation of nuclei used in heavy-ion collision experi- 
ments generates additional initial eccentricity (and consequently v 2 ) fluctuations on top of those 
existing in the collisions of spherical nuclei. 

1 To obtain non-zero contribution from b = Ofm collisions factor 2ir(b + 0.2)fm has been used. 
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IV Central collisions of deformed nuclei 



In the previous section we have shown that oblate and prolate deformations of nuclei generate 
significant fluctuations of the initial eccentricity at given impact parameter of collisions. 

What happens in the most central collisions ? To investigate this one has to calculate properly 
the centrality of collisions using number of participants N p3iTt and number of nucleon-nucleon 
collisions iV coll obtained from the optical Glauber simulation. We use two-component model [13] 
to calculate charged-particle multiplicity in the form: 



dN ch /drj 



[I 



AS 



x ) ■ n. 



part 



PP 



+ X 



(3) 



and quantity dN^/dr] is assumed to be the collision centrality measured experimentally. Results 
shown in Figure 3 have been obtained using n pp = 2.25 and x = 0.11 [IB]. Since the number N co \\ 
is significantly larger compared to N pait , approximately 40% of secondary particles originate from 
quantity iV co n and the rest of particle multiplicity is generated due to N pait . 
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Figure 3: Contour-plots of [ e ; dN^/drj ] values obtained with optical Glauber simulation 
for Au 197 +Au 197 collisions (left) and Ho 165 +Ho 165 collisions (right) at b = 0, 1, 2, . . . 12 fm. 



One observes (see also Figure 4) that most central collisions of oblate and prolate nuclei exhibit 
different behaviour. For prolate nuclei the highest multiplicity of secondary particles is predicted 
to happen when b = Ofm and spins of Ho 165 nuclei are parallel to the beam axis and to each 
other. In this case the number of nucleon-nucleon collisions is significantly higher compared to 
the case when spins of Ho 165 nuclei colliding at b = fm are orthogonal to the beams axis (and 
parallel to each other). Therefore we predict the elliptic flow values at very-high- multiplicity 
(VHM) collisions of prolate nuclei to decrease significantly (see Fig. 4). 

For oblate nuclei our optical Glauber model [7j predicts the opposite effect. The elliptic flow 
at most central (VHM) collisions should stay non-zero, and even slightly rise up ! This small 
increase of the eccentricity values in the highest multiplicity Au+Au collisions can be understood 
easily: For central collisions (6 = fm) the highest number of nucleon-nucleon collisions N co \\ is 
obtained when spins of two oblate nuclei are parallel to each other and orthogonal to the beam 
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axis. In this configuration eccentricity e of the interacting volume is non-zero (being proportional 
to deformation parameter j/^l of the oblate nuclei). Increasing the relative azimuthal angle of 
nuclear spins from — > 90 degrees in this configuration causes the eccentricity to vanish and 
number of n — n collisions N co \\ slightly decreases. 
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Figure 4: Contour-plots of [ e ; dN^/dr] ] values obtained for central Au 197 +Au 197 (left) and 
Ho 165 +Ho 165 collisions (right) at b = 0.0,0.2, . . . 3.2 fm. Error bars in bottom plots show the 
width of eccentricity fluctuations in the regions of dN^/drj indicated in contour-plots (top). 



This explains results shown in Fig.4 obtained for very high multiplicity Au+Au collisions 
simulated with impact parameters b = 0.0,0.2, . . .3.2fm. One should keep in mind that eccen- 
tricity values and fluctuations obtained by our optical Glauber model are subject to additional 
fluctuations originating from individual positions of interacting nucleons [H] which are taken into 
account only in Monte-Carlo Glauber model simulations [5]. Therefore, effects predicted here for 
Au+Au collisions may require high statistics and precise elliptic flow measurements at RHIC. 



V Comparison with RHIC data 

In heavy ion collisions experiments one can measure only the elliptic flow strength t> 2 (not the 
initial eccentricity e). However, hydrodynamical scenario [3] for the expansion of QCD matter 
created in these collisions predicts the elliptic flow to scale with initial eccentricity t>2 ~ k ■ e. 
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This allows one to compare relative strength of t>2 fluctuations a V2 / (t^) with the relative strength 
of eccentricity fluctuations a e / (e). 

In Figure 5 we compare results of Glauber Monte Carlo simulation (assuming Au 197 nucleus 
to be spherical) and experimental data obtained by PHOBOS collaboration [T3j with results of 
our optical Glauber simulation assuming Au 197 nucleus to be slightly deformed (/3 2 = —0.131). 

cf ( e ) / < e > [ NparTn 
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Figure 5: Relative fluctuations of eccentricity o(e)/(e) obtained in optical 
Glauber model (/?^ u = -0.13) in comparison with Glauber Monte Carlo 
simulation [T3] (assuming spherical Au 197 nucleus) and relative elliptic flow 
fluctuations a{v2)/{v2} as measured by PHOBOS collaboration [13] . 

We observe that elliptic flow fluctuations in the most central Au+Au collisions observed by 
PHOBOS collaboration can be almost fully accounted for assuming Au 197 nucleus to have oblate 
ground-state deformation \^\ ~ 0.13. Significantly lower strength of the elliptic flow fluctuations 
at smaller centralities predicted in our optical Glauber model for Au+Au collisions is a simple 
consequence of neglecting the eccentricity fluctuations originating from variations of the individual 
positions of interacting nucleons in the optical Glauber model. 

VI Conclusions 

Based on the results obtained with our simple optical Glauber model simulation we suggest 
to re-investigate consequeces of the possible ground-state nuclear deformation of Au 197 nucleus 
in Monte-Carlo Glauber model simulations. This would allow us to clarify the influence of the 
predicted oblate deformation of Au 197 nucleus on interpretations of relativistic Au+Au heavy ion 
collision experiments. We suggest to pay attention also to the most-central Cu+Cu collisions at 
RHIC since both stable Cu isotopes are predicted to be deformed [TU] . 

We have predicted that precise measurement of the elliptic flow and elliptic-flow fluctuations 
in most central nucleus-nucleus collisions may allow one to determine ground-state nuclear de- 
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formation of these nuclei. This should apply also to nuclei with zero spin (e.g. Si 28 , W 180 ) and 
unstable /exotic isotopes. 
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